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ABSTRACT 


In this paper we study the non-factorizable QED corrections to lU-pair-mediated (charged- 
current) four-fermion production in electron-positron collisions. A brief account of the ob¬ 
tained analytical results is given. They turn out to be different from the ones published in the 
literature. For the first time numerical results are presented, in particular the effects on the W 
line-shape. These effects are of the order of a per cent. Applying the same methods to ZZ- 
or ZH -mediated four-fermion production, the non-factorizable 0(a) corrections to the Z or H 
line-shape vanish. 
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1 Introduction 


With the start of LEP2, quantitative knowledge of the radiative corrections to the four-fermion 
production process e + e _ —> 4/ is needed [|T[. The full calculation of all these corrections 
will be extremely involved and at present one relies on approximations (p, such as leading- 
log initial-state radiation and running couplings [i;2j]. Another approach is to exploit the fact 
that in particular the corrections associated with the production of an intermediate IR-boson 
pair are important. This (charged-current) production mechanism dominates at LEP2 energies 
and determines the LEP2 sensitivity to the mass of the W boson and to the non-Abclian 
triple gauge-boson interactions. As such, one could approximate the complete set of radiative 
corrections by considering only the leading terms in an expansion around the W poles. The 
double-pole residues thus obtained could be viewed as a gauge-invariant definition of corrections 
to “VP-pair production”. The sub-leading terms in this expansion are generically suppressed 
by powers of Tw/Mw, with Mw and Tw denoting the mass and width of the W boson. The 
quality of this double-pole approximation degrades in the vicinity of the IF-pair production 
threshold, but a few Tjy above threshold it is already quite reliable [Q. It is conceivable that in 
the near future a combination of the above-mentioned approximations will result in sufficiently 
accurate theoretical predictions for four-fermion production processes. 

In the double-pole approximation the complete set of first-order radiative corrections to 
the charged-current four-fermion processes can be divided into so-called factorizable and non- 
factorizable corrections BB.i.e. corrections that manifestly contain two resonant W propaga¬ 
tors and those that do not. In view of gauge-invariance requirements, some care has to be taken 
with the precise definition of this split-up (see Sect. [2 . -1| ) . In the factorizable corrections one can 
distinguish between corrections to W- pair production and W decay. In this letter we address 
the question of the size of the usually neglected non-factorizable corrections. From the com¬ 
plete set of electroweak Feynman diagrams that contribute to the full O(a) correction, we will 
therefore only consider the non-factorizable ones, both for virtual corrections and real-photon 
bremsstrahlung. To be more precise, since we are only interested in the double-pole terms we 
are led to consider only non-factorizable QED diagrams in the soft-photon limit. Subsequently, 
the photon is treated inclusively, without imposing any limits on the photon phase space JiJ. 

This is the same approach as adopted by the authors of Ref. ||, who were the first to 
calculate non-factorizable W- pair corrections. For the present calculations, we have used two 
different methods. One is an extension of the treatment in the other is a modification of 
the standard methods, which involves a combination of the decomposition of multipoint scalar 
functions and the Feynman-parameter technique. The results obtained with our two methods 
are in complete mutual agreement. However, in contrast to 0 a clear separation between virtual 
and real photonic corrections has been made in both methods, which is essential to establish 
the cancellations of infrared and collinear divergences. This treatment reveals a significant 
difference between our results and those obtained by the authors of 0 . Our final results do not 
contain any logarithmic terms involving the final-state fermion masses, whereas in the results 
of 0 explicit logarithms of fermion-mass ratios occur (see discussion in Sect. 4.1 of Ref. |J). 
This difference can be traced back to the fact that although the fermion masses can formally 
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-Doi23 -Dq124 -E'01234 


Figure 1: Virtual diagrams contributing to the manifestly non-factorizable W- pair corrections. 
The scalar functions corresponding to these diagrams are denoted by -D 0123 , ^ 0124 , and -^ 01234 - 


be neglected in the absence of collinear divergences, they have to be introduced in intermediate 
results in order to regularize those divergences before dropping out from the final results. 

In the next section we briefly focus on the analytical results as obtained with the modified 
standard technique. A detailed account of our study and a discussion of both calculational 
methods will be published elsewhere. In the last section we present numerical implications 
of the non-factorizable corrections. Our calculations confirm that non-factorizable corrections 
vanish in the special case of initial-final state interference, thereby making non-factorizable 
radiative corrections independent from the W production angle, and in all cases when the 
integrations over both invariant masses of the virtual W bosons are performed ||. The practical 
consequence of the latter is that pure angular distributions are unaffected by non-factorizable 
O(a) corrections. So, the studies of non-Abelian triple gauge-boson couplings at LEP2 |J 
are not affected by these corrections. The completely new aspect that we have addressed in 
our analysis is the effect of non-factorizable corrections on invariant-mass distributions {W 
line-shapes). These distributions play a crucial role in extracting the PF-boson mass from 
the data through direct reconstruction of the Breit-Wigner resonances. The non-factorizable 
corrections to the line-shapes turn out to be the same for quark and lepton final states, provided 
the integrations over the decay angles have been performed. 


2 Non-factorizable corrections: analytical results 


In this section we present the results calculated in the modified standard technique. Details of 
the calculation and the alternative treatment, which is an extension of the method of [!|], will 
be published elsewhere. 
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2.1 Virtual corrections for purely leptonic final states 


As a first step we consider the manifestly non-factorizable corrections to the simplest class of 
charged-current four-fermion processes, involving a purely leptonic final state: 

e + (<?i)e - (g 2 ) —> PF + (pi) + W~{p 2 ) —> z/ £ (^)£ + (fci) + l'~{k 2 )vi>'{k' 2 ). (1) 

All external fermions are taken to be massless whenever possible. The relevant contributions 
consist of the final-final and intermediate-final state photonic interactions displayed in Fig. [[]. 
In principle also manifestly non-factorizable vertex corrections exist, which arise when the pho¬ 
ton in Fig. [I] does not originate from a PF-boson line but from the ^WW/ZWW vertex (hidden 
in the central blob). Those contributions can be shown to vanish, using power-counting argu¬ 
ments |T|. Also the manifestly non-factorizable initial-final state interference effects disappear 
in our approach. This happens upon adding virtual and real corrections. 


The double-pole contribution of the virtual corrections to the differential cross-section can 
be written in the form 


= 327tq: Re 


i(p 2 ■ A;i)DiAu 23 + D2D0124 + i(ki ■ k 2 )DiD 2 E 0 i 234: 


do 


Born? 


( 2 ) 


where D\ 2 = p\ 2 — M^ + iMw^w are the inverse (Breit-Wigner) PF-boson propagators. The 
functions ^ 0123 , An 24 , and £01234 are the scalar integrals corresponding to the diagrams shown 
in Fig. [Lj. with the integration measure defined as d 4 /c/(27r) 4 . The propagators occurring in 
these integrals are labelled according to: 0 = photon, 1 = PF + , 2 = IF - , 3 = A 1 ", and 4 = 
. Note that the factorization property exhibited in Eq. (§) is a direct consequence of the 
soft-photon approximation, which is inherent in our approach. 


To write down the analytical results we need to introduce some kinematic invariants: 

m i ,2 = ^ 1 , 2 ) s = (p 1 EP 2 ) 2 , S\ 22 ' = (Aq + k 2 + k ' 2 ) 2 , S 12 = (ki + k 2 ) 2 , (3) 

and some short-hand notations: 

Vo = ^ C = 1 - ^ - <0, x s = + *0, (3 = yjl-AMfa/s . (4) 

ffere ±*0 denotes an infinitesimal imaginary part, the sign of which is determined by causality. 


The scalar four-point function -D 0123 is infrared-finite, owing to the finite decay width of the 
PF boson. In the soft-photon limit it takes the form 


D, 


0123 


167 T [D 2 - C£>i] 


2 £i 2 


1 1 


2/o ’ C 


- £i 2 x s ] — - C\ 2 

Vo 


1 1 

Xs Vo 


+ £i 2 (x s ; ( j + £i 2 ( —; C ) + 


x x 


/ /If 2 \ 1 r 

ln (—r)+ 21n (C) Hvo)+HO 

V mi J \ L 


The function jC\ 2 (x]y) is the continued dilogarithm 

Ci 2 {x] y ) = Li 2 (l - xy) + ln(l - xy) [ln(xi/) - ln(x) - ln(t/)]. 


(5) 


( 6 ) 
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with Li 2 (x) the usual dilogarithm and x,y lying on the first Riemann sheet. The answer for 
the second four-point function, £ 0 i24, can be written in a similar way. 

The five-point scalar function, £'01234, can be evaluated by means of a decomposition into a 
sum of four-point functions 0 . In the double-pole approximation this decomposition reads 

w 1 £01234 = 2 A 4 £1234 + (w ■ Vi) £0234 + {w • U2) £0134 + {w • V3) £0124 + (w ■ V4 ) £0123, ( 7 ) 

with 


vin = v 2fl = +e afl p 1 p*k%kl, u 3m = - 

= +£ a p 1 nPiP2k'i, = £i< + £ 2 ^ 2 , A 4 = [£ a p lS Pil4kik s 2 ] 2 , (8) 

using the convention e 0123 = —£0123 = 1 - The labelling of the scalar functions (D, ]k i) is defined 
below Eq. (| 2 |). Note that the scalar four-point function £1234 is purely a consequence of the 
decomposition (| 7 ]). It does not involve the exchange of a photon and is therefore not affected 
by the soft-photon approximation. Since we are only interested in the double-pole residue, 
it should be calculated for on-shell W bosons. For the analytical expression, which is too 
involved to be presented here, we refer to the literature ||. The other new scalar four-point 
functions, £0134 and £0234, are infrared-divergent and should be calculated in the soft-photon 
approximation. Using a regulator mass A for the photon we can write 


£ 


0234 — 


167 T 2 Si 2 £ 


Li 2 (l + ^)-21nf^ 

S12 / V -£ 2 


In 


m\rri2 

— S 12 - *o 


7 r 

T 


i„ 2 (t!£ 

V m 2 


+ hr 


/ m 1 A 
\^Mw ) 


( 9 ) 


with a similar expression for £0134. 


2.2 Real-photon radiation for purely leptonic final states 


Only interferences of the real-photon diagrams can give contributions to the manifestly non- 
factorizable corrections. The relevant interferences can be read off from Fig. |l| by taking the 
exchanged photon to be on-shell. The infrared divergences contained in the virtual corrections 
will cancel against those present in the corresponding bremsstrahlung interferences. 

It should be noted that it is more complicated to obtain the five-point radiative interference 
correction. This is because the decomposition that we used in the case of the virtual five- 
point function cannot be carried over to the real-photon case. However, it is still possible to 
derive another decomposition using similar, but less straightforward arguments. Denoting the 
radiative analogues of the virtual scalar functions by a superscript ‘R’, we find 


w 


12 rriR 

A 


01234 


[w ■ rq) £q2 34 + [w ■ v 2 ) £0434 + {w ■ V 3 ) £q 424 
+ i w ‘ V4) ArL 23 + A 4 £^34. 


( 10 ) 
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The four-vectors w' and v[ are defined as before, but for real-photon emission. This is equivalent 
to the following substitutions: p\ — > —p\ , k\ —> —k\ and D 2 —■► D 2 . The radiation function 
D ^34 is an artefact of the decomposition fllCf) and does not involve the exchange of a photon. 
It can be obtained from .D1234 by the substitutions pi —> —pi and k\ —> —k±, resulting in the 
relation ImZ}^ 34 = I1T1.D1234. 


As will be explained in detail elsewhere, the radiative interferences can in fact be obtained 
from the virtual corrections by only considering the contribution from the photon pole to the 
complex k° integration and by making certain substitutions. The photon-pole part D^ kl of 
the scalar four-point function Diju. is obtained by subtracting the particle-pole part Df*‘ kl from 
^ k/jki ■ 

Dim = D ‘m ~ D <m- (U) 

The particle-pole parts are found to be 


7 -jpart 

-^0123 


7-xpart 

-^0234 


87 tMw D 2 - C-Di 


1 


1 


87 TS 12 D 2 


In 


Do 


iMlr 


ln(l - y 0 x s ) - ln(l - x s /Q 

_ _ hl (A) 


( 12 ) 

( 13 ) 


with similar expressions for -Dom and ^0134, respectively. The radiative interferences can be 
obtained from Eq. (Q) by adding a minus sign, by inserting the decomposition given in Eq. (^), 
and by substituting 


in the -D0123, -A134 terms: 
in the Aim, -A234 terms: 
in the -D1234 terms: 


-A)123, -Aim —i ► -Do 123, -^0134 
-Al24, D0234 > -Dqi24> -D0234 

D1234 —^ Df 234 


followed by 

D\ - 

+ —D* 

-^1, 

followed by 

d 2 - 

4. _n* 

^ J - y 2, 

followed by 

d 2 - 

+ —D* 

* 4 V 2 . 


2.3 Gauge-invariant definition of non-factorizable corrections 


The set of manifestly non-factorizable QED diagrams displayed in Fig. [I] is not gauge invariant. 
In order to achieve a gauge-invariant definition of the non-factorizable corrections, all (soft) 
photonic interactions between the positively (e + , IE + ,A) and negatively (e _ , W~ charged 
particles should be taken into account. Looking at Fig. |I} this is equivalent to the set of all up- 
down QED interferences. In the soft-photon, double-pole approximation only the “Coulomb” 
interaction between the off-shell W bosons survives as an extra contribution to the differential 
cross-section: 


dcr ?ivt(Pi\P2) = 327 TO: Re 


i{Pl ■ P 2 )C 0 12 


da 


Born* 


( 14 ) 


The scalar three-point function C012 is defined according to the notation of Sect. | 2 . 1 | . In our 
approximation it is artificially ultraviolet-divergent. Introducing an upper bound A for the 
energy of the photon, this scalar function reads 


0)12 


i 

16 tt 2 s /3 





2Li 2 ^l — —^ + - hr(i/ 0 ) 
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+ ln(z.) 


(15) 


In 


lD \ ) + In 
2M\yAJ 


iP‘2 \ 

2MwA/ 


— 2in In 


^)) 


The corresponding radiative interference can again be related to the virtual correction (B by 
adding a minus sign and by substituting C 012 —► Cq 12 followed by Di —> — D\. The photon-pole 
part Cqi 2 = C*oi 2 — Cq'i 2 can be derived from Eq. (|T% and 


^-ypart 

^012 


g^g|ln(l - x a ) + ln(l + a;,) - ln(l 


y 0 x s ) - In 


-iD 2 \ 
M W A] 


(16) 


If the virtual and real corrections are added, the dependence on the cut-off parameter A van¬ 
ishes. When we mention non-factorizable corrections in the following, we implicitly refer to the 
gauge-invariant sum of the manifestly non-factorizable corrections and the above-mentioned 
“Coulomb” contribution. 


2.4 Semi-leptonic and purely hadronic final states 


For the purely hadronic final states there are many more diagrams, as the photon can interact 
with all four final-state fermions. In order to make efficient use of the results presented in 
the previous subsections, we first introduce some short-hand notations based on the results 
for the purely leptonic (LL) final states. These short-hand notations involve the summation of 
virtual and real corrections to the differential cross-section. For instance, the virtual corrections 
originating from the first diagram of Fig. JT] can be combined with the corresponding real- 
photon correction into the contribution da^(kp, k[ \p 2 ). In a similar way virtual and real five- 
point corrections can be combined into da^l(kp, k[\k 2 ] k' 2 ). The gauge-restoring “Coulomb” 
contribution will be indicated by da c (pi\p 2 ). In terms of this notation the non-factorizable 
differential cross-section for purely leptonic final states becomes 

da LL (ki, k[\k 2 ] k 2 ) = da^kp k[\p 2 ) + da^l(k 2 , /4|pi) + da^l(h; k[\k 2 ; k 2 ) + da c (pi|p 2 ). (17) 


Analogously the non-factorizable differential cross-section for a purely hadronic final state ( HH ) 
can be written in the following way 


da H H(ki, k[\k 2 , k 2 ) = 3x3 ^ da^kp, k[\p 2 ) + | da^k'p, h\p 2 ) + ^da^l(k 2 , k 2 |pi) 

+ j d(T ( L 4 l(k' 2 ; k 2 \pi) + i i da^lih] k[\k 2 , k 2 ) + ^ ^ da^k'p, A; X |A; 2 ; k 2 ) 

2 da^lih; k[\k 2 ; k 2 ) + % ■ \ da^l(k[; ki\k' 2 , k 2 ) + da c (p]\p 2 ) 


3 3 


3 3 


(18) 


In order to keep the notation as uniform as possible, the momenta of the final-state quarks 
are defined along the lines of the purely leptonic case with ki (&') corresponding to down (up) 
type quarks. If one would like to take into account quark-mixing effects, it suffices to add the 
appropriate squared quark-mixing matrix elements (|W,-| 2 ) to the overall factor. Note that top 
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quarks do not contribute to the double-pole residues, since the on-shell decay W —> tb is not 
allowed. Therefore the approximation of massless final-state fermions is still justified. 

For a semi-leptonic final state (say HL ), when the W + decays hadronically and the PF~ 
leptonically, one can write 


da HL {kpk[\k 2 ;k 2 ) = 3 ^ da^kp, k[\p 2 ) + da^k'p, h\p 2 ) + da^l(k 2 ; k' 2 \pi) 
+ ^ da^lih] k[\k 2 , k' 2 ) + | K) + da c (p 1 \p 2 ) 


(19) 


Upon integration over the decay angles, the functions da z l and da^ L become symmetric 
under ki k[. As a result, the expressions (|l|) and (|1^) take on the form of (|T7| ) multiplied 
by the colour factors 9 and 3, respectively. These are precisely the colour factors that also 
arise in the Born cross-section. Therefore, after integration over the decay angles, the relative 
non-factorizable correction is the same for all final states. This universality property holds for 
all situations that exhibit the ki k[ symmetry. 


3 Numerical results 


In this section some numerical results will be presented. The quantity of interest is the relative 
non-factorizable correction 5 n f, defined as 


da 

d£ 


dUBorn 

d£ 


[1 + 


( 20 ) 


where £ represents some set of variables. Here we consider consecutively the distributions 
da/\dM\dM 2 d cos 6 1 ], da/[dMidM 2 \, da/dMi and da/dM aV} with Mi = , M av = |(Mi + 

M 2 ) and 6\ is the decay angle between k\ and p\ in the fab system. The results are shown in 
Fig. [2j for the angular distribution, and in Table [l] and Fig. [] for the invariant-mass distributions. 
The pure invariant-mass distributions play an important role in the extraction of the lU-boson 
mass from the data through direct reconstruction of the Breit-Wigner resonances. In this 
context especially the position of the maximum of these Breit-Wigner curves is of importance. 

All results in this section are presented for the following set of input parameters: 


M w = 80.22 GeV, T w = 2.08 GeV, M z = 91.187 GeV, r z = 2.49 GeV, 
a = 1/137.0359895, sin 2 6 W = 0.226074. 

From Fig. it is clear that corrections of a few per cent could arise for angular distributions. 
They should, however, vanish after integration over M\ and Ad 2 , as was mentioned before. The 
non-factorizable corrections <5 n /(Mi, M 2 ) to the double invariant-mass distribution are presented 
in Table [3]. From those results one can expect that 5 n f will be less steep for the Mi distribution 
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Figure 2: The relative non-factorizable correction M 2 , cos 9i) to the decay angular 

distribution da / [dMidM 2 d cos 61 ] for fixed values of the invariant masses M 12 [in GeV]. Centre- 
of-mass energy: y/s = 184 GeV. 


than for the M av distribution. This is confirmed by Fig. |3|. The corrections shown in Fig. ^ lead 
to a shift in the position of the maximum of the Breit-Wigner curves of the order of 1-2 MeV. 
These results have been obtained for the centre-of-mass energy y/s = 184 GeV. On the interval 
170-190 GeV the largest corrections are observed for 170 GeV, where the corrections are about 
a factor of two larger than those at 184 GeV. At 190 GeV the corrections are slightly smaller 
than those at 184 GeV. 


Ai 

a 2 


-1 

- 1/2 

- 1/4 

0 

1/4 

1/2 

1 

-1 

+ 0.84 

+ 0.65 

+ 0.51 

+ 0.37 

+ 0.21 

+ 0.06 

- 0.17 

- 1/2 

+ 0.65 

+ 0.52 

+ 0.39 

+ 0.24 

+ 0.07 

- 0.07 

- 0.26 

- 1/4 

+ 0.51 

+ 0.39 

+ 0.28 

+ 0.13 

- 0.02 

- 0.15 

- 0.31 

0 

+ 0.37 

+ 0.24 

+ 0.13 

+ 0.00 

- 0.13 

- 0.24 

- 0.37 

1/4 

+ 0.21 

+ 0.07 

- 0.02 

- 0.13 

- 0.24 

- 0.32 

- 0.43 

1/2 

+ 0.06 

- 0.07 

- 0.15 

- 0.24 

- 0.32 

- 0.39 

- 0.47 

1 

- 0.17 

- 0.26 

- 0.31 

- 0.37 

- 0.43 

- 0.47 

- 0.53 


Table 1: The relative non-factorizable correction 8 n f(Mi,M 2 ) [in %] to the double invariant- 
mass distribution da/[dMidM 2 ] for some particular values of M 12 . The invariant masses M l 2 
are specified in terms of their distance from Mw in units of Tw, he. Ai >2 = [Mi i2 — Mw\/Tw- 
Centre-of-mass energy: y/s = 184 GeV. 
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Figure 3: The relative non-factorizable correction S n f(M ) to the single invariant-mass distri¬ 
butions da/dM. Shown are the corrections to the distributions with respect to M\ and M av . 
Centre-of-mass energy: y/s = 184 GeV. 

4 Conclusions 


In this letter some analytical and numerical results are presented for non-factorizable corrections 
to VF-pair-mediated four-fermion production. In principle these corrections could be relevant 
for tests of triple gauge-boson couplings and for the determination of the IF-boson mass. For the 
latter the corrections are of O(a) and change the W line-shape by about 1%. For the former 
they vanish at the O(a) level. In view of the present experimental accuracy, the common 
practice of neglecting non-factorizable corrections is justified. 

One may wonder how non-factorizable corrections affect Z-pair-mediated and ZH -mediated 
four-fermion final states. In those cases, only five-point functions contribute, of which there 
are four contributions, as in Eq. ([T8|) . However, in contrast to Eq. fll8P , the charge factors are 
pair-wise opposite, such that integration over the decay angles leads to a vanishing result. Thus 
O(a) non-factorizable corrections to invariant-mass distributions in Z-pair-mediated or ZH- 
mediated four-fermion processes vanish. 


Acknowledgements: The authors are grateful to Dr. G.J. van Oldenborgh for useful 
discussions and for making some of his programs available to us. 
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